function prepare_graphs(file_ref,ref,output_directory,show_combined,show_income)

        %%%%%%%%%%%%%%%%%%%%%%%
        % LOAD DATA
        %%%%%%%%%%%%%%%%%%%%%%%

eval(['load ' output_directory 'data_' file_ref '.mat']);

        %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
        % FIGURES 
        %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

tt_vec=[0.25 0.5 0.75]';
pmin=0.2;
pmax=0.365;
qmin=[6.7 6.6 6.5;	
	  7.1 7.0 6.8;	
	  7.4 7.3 7.2];
qmax=qmin+0.4;

tt=2;
yy=2;

	prepare_graph_both     (tt,yy,p_interest,ghat_MLE_tau25_unconstr,ghat_MLE_tau25_unconstr_noBerk,ghat_MLE_tau25_constr,ghat_MLE_tau25_constr_noBerk,ghat_MLE_tau50_unconstr,ghat_MLE_tau50_unconstr_noBerk,ghat_MLE_tau50_constr,ghat_MLE_tau50_constr_noBerk,ghat_MLE_tau75_unconstr,ghat_MLE_tau75_unconstr_noBerk,ghat_MLE_tau75_constr,ghat_MLE_tau75_constr_noBerk,output_directory, ref,pmin,pmax,qmin,qmax,tt_vec,y_interest,y_interest_univ);

if show_income==1

	for yy=1:3

		prepare_graph_both_slim(tt,yy,p_interest,ghat_MLE_tau25_unconstr,ghat_MLE_tau25_unconstr_noBerk,ghat_MLE_tau25_constr,ghat_MLE_tau25_constr_noBerk,ghat_MLE_tau50_unconstr,ghat_MLE_tau50_unconstr_noBerk,ghat_MLE_tau50_constr,ghat_MLE_tau50_constr_noBerk,ghat_MLE_tau75_unconstr,ghat_MLE_tau75_unconstr_noBerk,ghat_MLE_tau75_constr,ghat_MLE_tau75_constr_noBerk,output_directory, ref,pmin,pmax,qmin,qmax,tt_vec,y_interest,y_interest_univ);

	end % yy

end

if show_combined==1	
	
	prepare_graph_combined(p_interest,ghat_MLE_tau25_unconstr,ghat_MLE_tau25_unconstr_noBerk,ghat_MLE_tau25_constr,ghat_MLE_tau25_constr_noBerk,ghat_MLE_tau50_unconstr,ghat_MLE_tau50_unconstr_noBerk,ghat_MLE_tau50_constr,ghat_MLE_tau50_constr_noBerk,ghat_MLE_tau75_unconstr,ghat_MLE_tau75_unconstr_noBerk,ghat_MLE_tau75_constr,ghat_MLE_tau75_constr_noBerk,output_directory, ref,pmin,pmax,y_interest,y_interest_univ);

end

	%

